Distributing fully optomechanical quantum correlations 



O 

(N 

oo 



Oh 



> 

00 

o 

^' 

o 



X 



L. Mazzola''^ and M. Patemostro^ 
' Turku Centre for Quantum Physics, Department of Physics and Astronomy, 
University of Turku, FI-20014 Turun yliopisto, Finland 
^ Centre for Theoretical Atomic, Molecular and Optical Physics, 
School of Mathematics and Physics, Queen's University Belfast, BT7 INN Belfast, United Kingdom 

We present a scheme to prepare quantum correlated states of two mechanical systems based on the pour- 
ing of pre-available all-optical entanglement into the state of two micro-mirrors belonging to remote and non- 
interacting optomechanical cavities. We show that, under realistic experimental conditions, the protocol allows 
for the preparation of a genuine quantum state of a composite mesoscopic system whose non-classical features 
extend beyond the occurrence of entanglement. We finally discuss a way to access such mechanical correlations. 
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I. INTRODUCTION 

Although there is no fundamental constraint to the achieve- 
ment of a quantum mechanical state of a macroscopic sys- 
tem, our daily experience is that we hardly observe any non- 
classical feature in intrinsically large or massive systems. A 
widely accepted justification for such difficulties relies on the 
paradigm of decoherence: it is not possible to perfectly iso- 
late a system from its surrounding environment. This de- 
termines an uncontrollable flux of information between such 
parties that acts to the detriment of any non-classicality the 
system was endowed with. Such effect appears to be severe 
for large-scale systems. Nonetheless, the spectacular level 
of sophistication reached in the manipulation and control of 
complex systems has now enabled the exploration of quan- 
tum mechanical effects on an increasingly large scale in vari- 
ous physical contexts. Such endeavors have produced impor- 
tant results: from the preparation of moderate-amplitude opti- 
cal Schrodinger cat states [1] to the engineering of multipar- 
tite entangled states in superconducting circuits [2], optically 
trapped atoms IJ] and Bose-Einstein condensates (j4l]. 

Recently nano/micro-scale vibrating structures have re- 
ceived considerable attention as potential platforms for the 
realization of ultra-sensitive measurements and as test-beds 
for probing the foundations of quantum mechanics. Passive 
and feedback-based cooling of micromechanical devices fsfc, 
strong optomechanical coupling [6-8] and, to some extent, 
the preparation of the ground state of a micro-mechanical 
mode jgl-lllll have been demonstrated in some outstanding ex- 
periments. An optimistic approach to the progress of this area 
of investigation makes the preparation of highly non-classical 
mechanical states a foreseeable possibility. This has trig- 
gered proposals to prepare optomechanical and fully mechan- 
ical entangled states |12-16], coherent macroscopic superpo- 
sitions ilTIl and squeezed states of mechanical modes OlSll 
achieved with either a 'gently' modulated driving or the use of 
non-classical light. A few steps towards hybridization through 
atomic -induced control of optomechanical systems have been 
performed [19]. 

In this paper we propose a scheme for the generation of a 
non-classical state of two spatially distant mechanical oscil- 
lators. The idea is to 'distribute' non-classical correlations 



from a pair of electromagnetic fields to two remote and non- 
interacting mechanical modes. This is done by feeding two in- 
dependent optomechanical cavities with a two-mode squeezed 
state of light [20]. We show how the light-matter interaction 
'transduces' the optical quantum correlations into fully me- 
chanical ones, robustly against the dissipative mechanisms af- 
fecting the setup, the mechanical modes' masses and the op- 
erating temperature. Besides the fundamentally interesting 
result of achieving entanglement between two massive ob- 
jects subjected to open-system dynamics, our proposal has 
also technological relevance. In fact, the scheme represents an 
experimentally viable strategy to distribute quantum correla- 
tions among spatially separated micro-mechanical devices. In 
this sense, our study embodies the continuous-variable coun- 
terpart of proposals for continuous-to-discrete-vaiiable entan- 
glement distribution [21], realized through the (for this task) 
largely unexplored mechanism of radiation pressure coupling. 
Moreover, our work goes beyond the proposal put forward in 
Ref. ifTill . where an adiabatic Hamiltonian model was con- 
sidered so as to 'transfer' the whole optical state onto the 
mechanical one. We achieve non-classicality in a much less 
stringent interaction regime, which makes our flexible proto- 
col different in principle and closer to the experimental reality. 

The basic idea of our scheme suggests also a way to indi- 
rectly infer the non-classicality of the mechanical modes. In 
fact, we show that the distributed fully mechanical entangle- 
ment can be mapped out onto the state of ancillary light modes 
and thus read through routine all-optical detection means. The 
feasibility of the proposal is verified using values for the rele- 
vant parameters in the problem boiTowed from a recent exper- 
iment [6], which makes our proposal a foreseeable possibility 
for building up a table-top network of optomechanical nodes 
connected by shared quantum correlations. 

The remainder of this paper is organized as follows. In 
Sec. HI] we introduce the system at hand and describe our pro- 
posed protocol for entanglement distribution. Sec. |III]is de- 
voted to the quantification of the non-classical correlations in 
the purely mechanical system. Such an analysis will be per- 
formed by determining both entanglement and more general 
quantum correlations. Sec. llVldesciibed our proposal for en- 
tanglement inference, which is based on the use of a pair of 
ancillary optical modes. Finally, in Sec.|V]we draw our con- 
clusions. 
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FIG. 1: (Color online) Sketch of the thought experiment: two 
optomechanical cavities are pumped by classical laser fields (not 
shown) and two-mode squeezed light, generated by feeding a non- 
linear crystal. We show the symbols for polarization beam split- 
ters (PBSs) and quarter wave-plates (QWPs) that are used to direct 
the polarized input squeezed light (output field) to the cavity (the 
homodyne-detection stage). 



II. THE MODEL AND PROTOCOL 

We consider two independent and non-interacting optome- 
chanical systems, each consisting of a Fabry-Perot cavity 
composed of a heavy fixed input mirror and a lighter movable 
end one. The movable mirrors are modelled as single-mode 
quantum harmonic oscillators coupled via radiation-pressure 
to the corresponding intracavity field (see Fig. 1). Each cavity 
has length L, and frequency wj-, and is pumped by a laser field 
with frequency ljl. In addition, the two cavities are jointly 
pumped by a two-mode light field (each mode has frequency 
ojs ), prepared in an entangled state. The Hamiltonian of the 
two optomechanical devices, in a frame rotating at the fre- 
quency of the lasers, reads 
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where (p,) is the position (momentum) operator of the /* 
mirror, c, (c] ) is the annihilation (creation) operator of the /* 
cavity field (whose photon-number operator and energy decay 
rate are n, and /c,, respectively), 6'=a)'^-a>L is the pump-cavity 
detuning. Moreover, wj,, and m, are frequency and mass of 
the /* mechanical oscillator and;^, = wJ-,/L, is the coupling 
rate with the corresponding cavity field. Finally, we have £, - 



^^^^ ^1.2 the power of the pumping fields. 



The system is affected by cavity losses and decoherence 
induced by the thermal Brownian motion of the mechanical 
oscillators, which are damped at a rate yj^. The model em- 
bodied by Eq. ([TJ encompasses a dynamics of non-trivial so- 
lution due to the explicitly nonlinear nature of the radiation- 
pressure term. A considerable simplification, though, comes 
from having intense fields feeding the cavities, thus allowing 
us to linearize the fields' and mirrors' operators around their 
respective steady-state. The equations of motion can then be 
cast into the compact form 120,1 



where fJ-{5Qi, dPj, (5x,, Syi) is the ordered vector of the fluc- 
tuations of the dimensionless quadrature operators Qi = 
qi -ijmicij'fjjh and P,- - pj yjtmiiLo'^ for mechanical mode / and 
6xi = {Sc:] + 5c,)/ V2, 6yi = i{6c] - 6ci)/ V2 for the corre- 
sponding cavity fields. Each 4x4 kernel matrix K, reads 
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with gi=Xi V^/(2'W|Wj„) being an effective coupling rate, 
c'^-&i I {ki + iAi) the amplitude of cavity field / and 
Ai-<J^-(jJi-Xq's '^he corresponding effective cavity-laser de- 



tuning. Finally q'^ 



.MM!! 



is the steady state displacement 



of the mechanical mode /. The last term in Eq. (|2]l is the 
vector of input noise nj=(0,|,-, y/lKidx'.^, y/lKidy'jJ, where 

I, is the zero-mean Langevin force operator accounting for 
the Brownian motion affecting the mechanical mode /. For 
large mechanical quality factors, such term is correlated 
as {l(t)l(t'))=(2yi„kBTi)6(t-t')/J,„, with kg the Boltzmann 
constant and T, the temperature of the /* mechanical bath, 
while 6^.„HStl+6c'.^/ V2 and 6y'.„^i(6c'l - 6c'.)l ^ are the 
quadratures of the input noise to a cavity. 

Clearly, the dynamics encompassed by Eq. ^ would give 
rise to two independent evolutions, one for each optomechan- 
ical device. This seemingly prevents the establishment of 
any quantum correlation between the two mechanical modes, 
unless a pre-available resource is used. In fact, as proven 
in Refs. lf2ll I22I1 . the sharing of an off-line prepared en- 
tangled resource can be used in order to set entanglement 
between remote and non-interacting subsystems, subjected 
to bilocal interactions with the components of such entan- 
gled aid. We thus consider the two cavities pumped by 
non-classically correlated light as embodied by a two-mode 
squeezed vacuum state i23[l . Each mode of such pump has 



frequency (jj^-^L+^m (where we take equal mechanical fre- 

1 2 

quencies w,„ = w,„' ) so that the corresponding input-noise 
correlations C; = {{8c^'^ 5c'. ), {5c'. 5c'.'^'), {6c'. 5c': ), <5c^'"'"5c*"'")) 

J in m'' ^ m m'' ^ in m'' ^ in in ' ^ 

[for ii^k-1,2] can be expressed as Cj - R(t,t')6(t - f) 
with R(t,t')=(N,N + l,e-'""'<'-'')M,e''^"'<'-'''M*) a vector de- 
termined by the parameters sinh^ r and M- sinh r cosh re"^ 
that characterize the squeezed state (r and are the modulus 
and phase of the squeezing parameter). From now on, without 
affecting the generality of our study, we take = 0. 

Although the analytic expressions of the elements of f, can 
be found straightforwardly, they are too lengthy to be reported 
here. However, they can be handly used so as to calculate the 
correlation functions of the relevant operators of the system as 



JdojJ t/Oe-'*'^+"^'y„/j(w,Q) (4) 



5,f,- = K,f; + fi;, (/=1,2) 



(2) 



with Vap(cj,Q.) = {{g,(cj),gfim}/2 (a,jS=l,..,8) die 
frequency-domain correlation function between elements a 
and y6 of g-i5Qi, 6Pi , 5Q2, 6P2, 5xi , 6yi , 6x2, 5y2)- The matrix 
V(f) with elements defined as in Eq. (|4| embodies the time- 
dependent covariance matrix of the whole system. This is an 
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FIG. 2: (Color online) Logaritmic negativity against the effective 
detuning A and the dimensionless squeezing parameter r for T = 
2 mK. The remaining set of parameters is borrowed from the recent 
experiment reported in Ref. |@]. 



important remark as the enforced linearity of the model treated 
here guarantees that the dynamical map arising from the solu- 
tion of Eq. dU preserves the Gaussian nature of any input state. 
We consider a Gaussian input noise (the two-mode squeezed 
vacuum state) and a thermal state of the mechanical modes 
determined by the quantum Brownian motion. As Gaussian 
states are fully determined by their associated covariance ma- 
trix, the convenience of such description is apparent. 



III. NON-CLASSICAL CORRELATIONS 

We are now in a position to assess the non-classical cor- 
relations settled between the two mechanical modes. We do 
this using a twofold approach: first we quantify the entan- 
glement shared by mechanical modes 1 and 2, demonstrating 
the effectiveness of the process put forward here. Second, we 
demonstrate that the parameter-domain where purely mechan- 
ical non-classical correlations can be observed extends far be- 
yond the one where entanglement is non-zero. Let us start 
studying the degree of entanglement between the mechanical 
modes. Clearly, while the achievement of a non-zero degree 
of entanglement would show the plausible nature of our pro- 
posal as a way to distribute mechanical entanglement across 
a network of optomechanical nodes, its effectiveness will be 
evaluated in terms of the ratio between the set mechanical en- 
tanglement and the amount consumed as a resource. We use 
the logarithmic negativity as a measure of entanglement in 
the state of our system. Given the covariance matrix Vm of 
the mechanical system (obtained from V by taking the first 
four rows and columns), the logarithmic negativity is found 
as £■= max[0, - ln2(v_)] with v_ the smallest eigenvalue of 
the symplectic spectrum of = PV,„P with V-t®cr~ {cri is 
the ; — x,y,z Pauli matrix) [24]. 

As we aim at proposing a viable strategy for the settlement 
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FIG. 3: (Color online) (a) Logaritmic negativity against the dimen- 
sionless squeezing parameter r and the temperature of the mechan- 
ical bath T (in Kelvin) for A = lo„,. (b) Gaussian quantum discord 
against squeezing and temperature for A = oj,,,. Non-classicality of 
the mechanical modes is found up to T = 0. 12 K. 



of all-mechanical entanglement, it is important to test our pro- 
tocol against cuiTent experimental capabilities. Therefore, 
in our simulations we have used parameters taken directly 
from a very recent experiment [6|. In particular, we have 
taken identical mechanical modes with w„,/27r = 947 KHz, 

1 2 

mi. 2 = 145 ng and I2n - 140 Hz and cavities hav- 
ing length Li,2 - 25 mm, wavelength 1064 nm, decay rate 
A-/27r = 215 KHz and pumped by laser fields of 11 mW 
power Ii25i] . The free variables in our simulations are the ex- 
perimentally tunable detuning A; 2=A, the squeezing param- 
eter r and the initial temperature of the mechanical mirrors. 
Fig.|2]shows the logarithmic negativity for T - 2 mK and re- 
veals that the mechanical entanglement is maximum close to 
the frequency Um and non-zero in a region of A that is wide 
roughly k. A non-trivial dependence on r is also observed. 
As the squeezing increases, the region of values of A where 
E naiTows. At the same time, E first slowly increases 
with r, reaching a maximum at r ^ 1 and quickly decreas- 
ing for larger values {E disappears completely for r > 1.6). 
This behavior reminds of the analogous features observed in 
the transfer of entanglement to two remote qubits [21]. How- 
ever, while in that case the responsible for such features was 
the mismatch in the Hilbert-space dimensions, here the expla- 
nation is that with increasing squeezing also the thermal noise 
entering each cavity is enhanced [26]. We emphasize that two- 
mode squeezed states with squeezing parameters of the order 
of unity are well within experimental reach, as demonstrated, 
for example, in Ref. [ZT^]. 

We now investigate the behaviour of mechanical non- 
classical correlations beyond entanglement, as assessed by 
Gaussian quantum discord Ii28..1 . Quantum discord is defined 
from the discrepancy that two classically equivalent expres- 
sions of mutual information attain when extended to the quan- 
tum domain [29]. It is connected to the information on a 
system that cannot be retrieved by performing measurements 
only on a part of it. For Gaussian continuous-variable sys- 
tems, the evaluation of discord relies on an optimization pro- 
cedure over the set of any possible single-mode Gaussian 
measurements. For a two-mode Gaussian state having covari- 
7 
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trices accounting for the local variances of mode 1 (2) and 
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FIG. 4: (Color online) Entanglement input-output relation showing 
the mapping of mechanical entanglement onto optical one shared by 
the auxiliary modes. The quantities on the horizontal and vertical 
axes are dimensionless. 



the inter-mode correlations, the Gaussian quantum discord 
is £) = /( VAl) - f(fi-) - f{^i^) + inf^„ /( Vd^). Here, 
/(x) = (^)log[^] - (^)log[^], A2 = deta2, A<± are 
the symplectic eigenvalues of cr, e - a.^ - 7(a2 + co) '7^ 
is the Schur complement of ai, and ctq is a single-mode ro- 
tated squeezed state. Fig. |4]compares the behavior of E (part 
(a)) and D (part (b)) for the two mechanical modes when 
studied against T and r for A = u,,,. The two indicators of 
non-classicality reveal rather different trends when analyzed 
against the squeezing parameter. At low temperatures, the 
maximum of O is attained for r ^ 1, while for r>10mK 
(when entanglement disappears), D grows with the squeez- 
ing. Remarkably non-classical correlations still persists and 
are significantly non zero for temperature higher than O.IK, 
thus demonstrating that the model addressed in our proposal 
is able to set robust and general non-classical correlations in 
quite an ample region of the working parameters. 



IV. ENTANGLEMENT INFERENCE 

Here we put forward a proposal to reveal the entanglement 
set between the mechanical modes. Our scheme relies on the 
use of two auxiliary light modes, each prepared in a coher- 
ent state and interacting with the respective optomechanical 
cavity. The mechanical modes are assumed to be prepared 
in an entangled state of the sort addressed above. The con- 
cept at the basis of our inference scheme is similar to the one 
underlying the entanglement-setting one: the optomechanical 
interaction acts as a transducer, this time changing entangle- 
ment from mechanical to fully optical. As the ancillary modes 
are initially uncorrected and do not interact directly, the cre- 
ation of optical entanglement can only be a signature of the 
pre-existence of mechanical one. Differently from [IS], our 
scheme does not require an adiabatic light-mirror interaction 
nor any weak coupling approximation between mechanical 



and photonic ancillary modes. 

We take a fully dynamical approach based on the solution 
of the equation satisfied by the covariance matrix V of the me- 
chanical modes and the auxiliary fields, which can be shown 
to read 5,V = LV + \\J + D. Here L is the kernel of the 
Langevin equations of each ancilla-mirror system (L= ©^^j K, 
if we take the same parameters of the intracavity modes for the 
auxiliary ones), and D - Diag[0, y,„(2n + l),/c, /c, 0,7m(2n + 
\),k,k\ is the noise matrix with h - {exp{fi(L>„/kBT} - 1)"'. 
Figure [3] shows the conversion of input mechanical entan- 
glement in output optical one at A = cjjn^, T - 2 mK and 
r e [0, 2]. Remarkably, the output optical entanglement turns 
out to be linearly dependent on the input mechanical one. As 
the former can be determined by means of homodyne mea- 
surements on the ancillary modes, this suggests the possibil- 
ity to quantify mechanical entanglement by subjecting the re- 
constructed all-optical one to a post-processing procedure that 
rescales it according to the linear behavior highlighted above. 



V. CONCLUSIONS 

We have introduced a protocol for the engineering of spa- 
tially distributed mechanical non-classicality through a sim- 
ple bi-local radiation-pressure mechanism. By 'consuming' 
a pre-available resource of entanglement, we have achieved a 
highly non-classical state of two mechanical oscillators, suit- 
able for the distribution of entanglement and general quantum 
coiTelation in a network of continuous-variable systems. We 
have examined both the stationary and fully time-dependent 
dynamics of the system, finding robustness of the enforced 
quantumness in a wide range of operating parameters. We 
have also discussed a rather promising scheme for the infer- 
ence of the distributed fully mechanical quantum correlations 
based on a simple all-optical detection scheme and straight- 
forward post-processing of data. Although only hinted in our 
work, it will be very interesting to explore deeply the potential 
of our scheme as a mean to merge optomechanics and the the- 
ory of quantum memories |30] towards the construction of a 
quantum network whose nodes are embodied by high-quality 
mechanical system ll3ll[32|] . 
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